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Abstract. Three-body resonances in atomic systems are calculated as complex- 
energy solutions of Faddeev-type integral equations. The homogeneous Faddeev- 
Merkuriev integral equations are solved by approximating the potential terms in a 
Coulomb- Sturmian basis. The Coulomb-Sturmian matrix elements of the three-body 
Coulomb Green's operator has been calculated as a contour integral of two-body 
Coulomb Green's matrices. This approximation casts the integral equation into a 
matrix equation and the complex energies are located as the complex zeros of the 
Frcdholm determinant. We calculated resonances of the e — Ps system at higher 
energies and for total angular momentum L — 1 with natural and unnatural parity. 
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1. Introduction 

The wave function of a three-particle system is very complicated. It may have several 
different kinds of asymptotic behavior reflecting the possible asymptotic fragmentations. 
It is very hard to impose all the asymptotic conditions on a single wave function. The 
Faddeev approach is a simplification: the wave function is split into components such 
that each component describes only one kind of asymptotic fragmentation [1J. Then 
only one kind of asymptotic behavior should be imposed on each component. The 
components satisfy a set of coupled equations, the Faddeev equations. 

If we want to apply this idea to systems with Coulomb potentials, we may run 
into difficulties. The Coulomb potential is a long range potential, thus the motion in a 
Coulomb field never becomes a free motion, even at asymptotic distances. Consequently 
the separation of the wave function along different asymptotic properties does not really 
work. If we just plug the Coulomb potential into the original Faddeev equations, the 
equations become singular. The usual asymptotic analysis fails to provide the boundary 
condition. In integral equation form, the kernel of the equations fails to be compact 
and we cannot approximate them by finite-rank terms. Merkuriev proposed [U [2] a 
modification of the Faddeev procedure which led to integral equations with compact 
kernels and differential equations with known boundary conditions. 
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Resonances are related to the outgoing-wave solutions of the Schrddinger equation 
at complex energies E = E r — iT/2. Here E r is the resonance energy and T is 
the resonance width, which is related to the lifetime of the decaying state. In an 
integral equation formalism, the resonances are the solutions of the homogeneous integral 
equations on the unphysical sheet, close to real energies. 

A few years ago, a method for solving Faddeev-type integral equations for scattering 
problems [3] was adopted to calculate resonances [3J. The homogeneous version of the 
Faddeev integral equations was solved at complex energies on the unphysical sheet. The 
method entails expanding the potentials terms in the integral equations on a Coulomb- 
Sturmian basis. This transforms the integral equations to a matrix equation. S"-wave 
resonances of the e — Ps system have been calculated and good agreement with the 
results of methods based on complex rotation of the coordinates were found E] . 

An accumulation of resonance poles around thresholds has been reported in Refs. 
[7], however, we are not investigating these threshold resonances here. They are either 
too close to the threshold, too close to each other, or too broad as we move away from 
the threshold. In any case, their experimental verification does not seem to be likely in 
the near future. 

In this paper, we will report some new developments of this method. In Sec. 2, we 
will outline the Faddeev-Merkuriev approach to the three-body Coulomb problem. In 
Sec. 3, we detail the Coulomb-Sturmian separable expansion approach. We introduce 
a new contour integral for the three-body Green's operator which makes the analytic 
continuation to the resonance region easier. In Sec. 4 we present resonances of the e — Ps 
three-body system up to the fifth threshold with total angular momentum L = 1. We 
provide results for both natural and unnatural parity states. 

2. Faddeev-Merkuriev integral equations 

The Hamiltonian of an atomic three-body system is given by 



where H° is the three-body kinetic energy operator and denotes the Coulomb 
interaction of each subsystem a = 1,2,3. Throughout, we use the usual configuration- 
space Jacobi coordinates x a and y a , where x a is the distance between the pair 7) 
and y a is the distance between the center of mass of the pair ((3,j) and the particle 
a. Thus, the potential v%, the interaction of the pair ((3,j), appears as v^(x a ). In 
an atomic three-body system, two particles always have the same sign of charge. So, 
without loss of generality, we can assume that they are particles 1 and 2, and therefore 
v$ is a repulsive Coulomb potential. 

The Hamiltonian ([T| is defined in the three-body Hilbert space. Therefore, the 
two-body potential operators are formally embedded in the three-body Hilbert space, 



H = H +Vi +v% + v\ 



(1) 




(2) 
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where l Va is a unit operator in the two-body Hilbert space associated with the y a 
coordinate. 

The role of a Coulomb potential in a three-body system is twofold. The Coulomb 
potential is a long range potential but it also possesses some features of a short- 
range potential. It strongly correlates the particles and may even support two-body 
bound states. These two properties are contradictory and require different treatment. 
Merkuriev proposed a separation of the three-body configuration space into different 
asymptotic regions [2]. The two-body asymptotic region Q a is defined as a part of the 
three-body configuration space where the conditions 

(\x a \/x y < \y a \/y , (3) 

with parameters Xq > 0, yo > and v > 2 are satisfied. It was shown that in Q a the 
short-range character of the Coulomb potential prevails, while in the complementary 
region the long-range character of the Coulomb potential becomes dominant. Thus, it 
seems to be a good idea to split the Coulomb potential in the three-body configuration 
space into short-range and long-range terms 

= (4) 

where the superscripts s and I indicate the short- and long-range attributes, respectively. 
The splitting is carried out with the help of a splitting function ( a , 

Va\x a ,y a ) = v^(x a )(^ a (x a ,y a ), (5) 
v®(x a ,y a ) = v°(x a ) [1 -( a (x a ,y a )] . (6) 

The function ( a vanishes asymptotically within the three-body sector, where x a ~ y a — > 
oo, and approaches 1 in the two-body asymptotic region fl a , where x a « y a — > oo. As 
a result, in the three-body sector, v$ vanishes and v$ approaches v%. In practice, the 
functional form 

Ca(x a , y a ) = 2/ {1 + exp [(x a /x ) v / (1 + y a /y )]} (7) 

is used. Typical shapes for and are shown in Figures [l] and [2J respectively. 

In the Hamiltonian ([IJ the Coulomb potential v% is repulsive and does not support 
bound states. Consequently, there are no two-body channels associated with this 
fragmentation and the entire v% can be considered as long-range potential. Then the 
long-range Hamiltonian is defined as 

H® = H° + v®+v!!>+v$, (8) 

and the three-body Hamiltonian takes the form 

H = H®+v[ a) +vi'\ (9) 

This Hamiltonian looks like an ordinary three-body Hamiltonian with two short range 
interactions. 

To determine the bound and resonant states, we have to solve the Schrodinger 
equation 

H\9)=E\V) (10) 
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Figure 1. Short range potential t>W for attractive Coulomb potential. The parameters 
are Z = —1, x$ — 30, yo — 35 and v = 2.1 (in atomic units). 




Figure 2. Long range potential for attractive Coulomb potential. The parameters 
are as in Fig. 1. 
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for real and complex E eigenvalues, respectively. In the Faddeev approach the Faddeev 
components are defined by 

\1> a ) = (E-H®)- 1 vU\*), (11) 

where a = 1, 2. This involves a splitting of the wave function into two components 

I*) = |V>i> + |V>2>- (12) 
Then for the Faddeev components, we have the set of equations, the Faddeev equations, 

(E-H?)\^)=v{ s) \i; 2 ) (13) 

(E-H?M 2 )=v { 2 s) \^), (14) 

where 

ff(Q = #(0 + v W. (15) 



By adding these two equations and taking into account Eq. (12) we recover the original 
Schrodinger equation. So, the Faddeev procedure is no more and no less than a method 
of solving the Schrodinger equation. We can cast these differential equations into an 
integral equation form 

\i> 1 ) = G<t ) (E)vi a) \i> 2 ) (16) 
|^> = G?(£)4 s) |V>i>, (17) 

where G%\e) = {E-H i J ) y 1 . 

Before going further, we should examine the spectral properties of the Hamiltonian 

iZf = ff(0 + v { s) =H° + v? + v® + vl (18) 

It is obvious that it supports infinitely many two-body channels associated with the 
bound states of the attractive Coulomb potential vf ■ Potential v% is repulsive, therefore 
does not support bound states and there are no two-body channels associated with 
fragmentation 3. The three-body potential v% is attractive. It is a valley along a 
parabola-like curve which becomes shallower and shallower, and finally disappears as y 2 
goes to infinity (see Fig. 2). Thus, v 2 (^2,2/2) does not support two-body bound states 
either in the subsystem x 2 if y 2 00. Consequently, there are no two-body channels 
associated with fragmentation 2. Therefore, the asymptotic Hamiltonian H± has two- 
body channels only in the fragmentation where particle 1 is at infinity and particles 2 
and 3 form bound states. If either particle 2 or 3 is at infinity, no bound states are 
allowed in the respective subsystem. The corresponding Gf Green's operator, acting 



on the v[ \ip 2 ) term in (16), will generate only those type of two-body channels in \tpi) 
where particle 1 is at infinity and particles 2 and 3 form bound states. A similar analysis 
is valid also for l^)- So, the Merkuriev procedure results in a separation of the three- 
body wave function into components in such a way that each component has only one 
type of two-body channel. This is the main advantage of the original Faddeev equations 
and, as the above analysis shows, this property remains valid also for attractive Coulomb 
potentials. 
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The long-range part of the Coulomb potential, u« , does not support two-body 
channels. It may, however, support bound states, i.e. H"> may have three-body bound 
states. This can lead to the appearance of spurious solutions of the Faddeev-Merkuriev 
equations. If H"> has a bound state, then (z — if^) -1 is singular at this energy. 



Consequently, in Eq. (11), applying this singular operator on a vanishing may 
produce a non-vanishing \ip a ). So, we may find a solution where neither nor \ip 2 ) 
vanish, but = l^i) + IV^) vanishes. These states would be non-trivial solutions of the 
Faddeev-Merkuriev equations, but would be trivial solutions of the original Schrodinger 
equation. These states are spurious, or ghost, solutions. A way to eliminate them is to 
ensure that in the energy range of physical interest does not have bound states. We 
can achieve this by choosing the parameters Xq and yo accordingly. For resonances at 
higher energies we should take a bigger xq, thus pushing the unwanted bound states of 
out of the spectrum of physical interest. It is also obvious from this analysis that 
the spurious solutions are sensitive to the choice of xq and yo, while the true resonances 
are not. By varying the parameters x and y , one can single out the possible spurious 
solutions. 

A very nice advantage of the Faddeev equations is that the identity of particles 
simplifies the equations. If particles 1 and 2 are identical particles, the Faddeev 
components and \ip2), m their own natural Jacobi coordinates, must have the same 
functional forms 

(xm\ih) = (salted)- (19) 
On the other hand, by interchanging particles 1 and 2, we have 

VuM=pM, (20) 

where p — ±1 . Building this information into the formalism, we arrive at the integral 
equation 

= G?v{ s) pV 12 \^), (21) 

which by itself determines We notice that so far no approximation has been made, 
and even though this integral equation has only one component, it gives a full account 
of the asymptotic and symmetry properties of the system. 

3. Separable expansion solution of the Faddeev equations 

3.1. Coulomb- Sturmian basis 

The Coulomb-Sturmian (CS) functions [Sj are the solutions of the Sturm-Liouville 
problem of the Coulomb Hamiltonian 

d 2 1(1 + 1) 2b(n + l + l) A , , , 

+ ^—T^ L + b 2 )(r\nl;b)=0, (22) 



dr 2 r 
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where 6 is a parameter, n is the radial quantum number and I is the angular momentum. 
In configuration space, the CS functions are given by 



(r\nl; b) = ^ {n + ^ + ^ exp(-6r)(26r)' +1 4 2 ' +1 )(26r) , (23) 

where L denotes the Laguerre polynomials. By defining the functions (r\nl;b) = 
(r\nl;b)/r, the orthogonality and completeness relations take the forms 

(pfl^b\nl\ b) = (n'l; b\rrifb) = 5 nn , (24) 

and 

N N 

1= lim y^\nl;b)(nl;b\ = lim V \nl; b)(nl; b\ . (25) 

N^oo ' ^ N—>oo ' ^ 

n=0 n=0 

Since the three-body Hilbert space is a direct product of two-body Hilbert spaces, an 
appropriate basis is the bipolar basis, which can be defined as the angular-momentum- 
coupled direct product of the two-body bases, 

\nul\;b x b y ) a = \nl;b x ) a (g>\uX;by) a , (n, v = 0, 1, 2, . . .), (26) 

where \nl; b x ) a and |z/A; b y ) a are associated with the coordinates x a and y a , respectively. 
With this basis, the completeness relation takes the form (with angular momentum 
summation implicitly included) 

N 

1= lim \ \nisl\;b x b y ) a a (ni/l\;b x b y \ = lim 1^, (27) 

n,u=0 

where (xy\ni / lX]b x by) = (xylnulX; b x b y )/(xy). 
3.2. Separable approximation 

We may introduce an unit operator into the Faddeev equation 



|Vi> = lim G«\E)l?v{ s) l»\ilj 2 ) (28) 



\fa)= lim G^)WWi>- (29) 
This identity becomes an approximation if we keep N finite, which is the equivalent of 

(s) 

approximating v a in the three-body Hilbert space by a separable form 

_ l im -iN AsUN 1 JV (s)-,N 

N— too 
N 

\nul\, b x b y ) a v^ p(n'v'V\'; b x b y \, (30) 

n,v,n' ,v'=0 

where v^l = a (nvl\;b x b y \va S> \n'i''l'\';b x by}{3. These matrix elements can be evaluated 
numerically by using the transformation of the Jacobi coordinates [9j . The completeness 
of the CS basis guarantees the convergence of the expansion with increasing N and 
angular momentum channels. 
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Now, by applying the bra (n"i/"l"X";b x b y \ from the left, the solution of the 
homogeneous Faddeev-Merkuriev equation turns into the solution of a matrix equation 
for the component vector 



where 



G® = a (nulX;b x b y \G^\n'u'l'X';b x b y ) a . 
These equations can be transformed into the matrix form 



1 
1 



G?(E) 







Ul2 





0. 



(31) 
(32) 

(33) 
(34) 



G_2 (E) / \ i_ 2J _ 
which exhibits a homogeneous algebraic equation for the Faddeev components. This 
homogeneous algebraic equation is solvable if and only if 



D(E) = det 



g±\e) 
o gf{E) 



)■'- 




V.12 \ 






o J 



0. 



(35) 



where D(E) is the Fredholm determinant. The real-energy solutions provide us with 
the bound states, while the complex-energy ones give the resonant states. 

3.3. Calculation of G® 

Unfortunately, the Green's operator is not known. It is related to the asymptotic 
Hamiltonian H a l \ which is still a complicated three-body Coulomb Hamiltonian. 
However, Ha has only one type of two-body asymptotic channels where particle a 
is at infinity. This asymptotic Hamiltonian is denoted by H® s . If a three-body system 
has only one type of asymptotic channel, then a single Lippmann-Schwinger equation 
provides a unique solution: 

G${z) = G a a s {z) + G a a s (z)V: s G«\z), (36) 

where G\ s is an asymptotic channel Green's operator, G^ s (z) — (z — H^)' 1 , and 
V£ s = — H^ s . Merkuriev constructed G^ s in the different asymptotic regions 
of the three-body configuration space and showed that the kernel of this Lippmann- 
Schwinger equation is completely continuous (compact) [TJ [2] . Therefore, V£ s can also 
be approximated by a separable form 

iN-iras-tN ^ -tN-iras-tN 



lim i"v as l" « l n V as r 

„ a a a a a a 

N 



\nvl\\ b x b y ) a V^ a (n'u'l'X'; b x b y \, 



(37) 



n,u,n'u'=0 



where V^ = a {nulX] b x b y \V^ s \n'u'l'X'; b x b y ) a . The solution of Eqs. (36) can be expressed 
formally as 

(G^r 1 = (git 1 - n s , (38) 
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where 

G a a s = a (nvl\) b x b y \G a :\n'u'l'\'- b x b y ) a , (39) 

V a a s = a (nul\; b x b y \V: s \n'u'l'X'; b x b y ) a . (40) 



The matrix elements (|39j) and (40) have to be calculated between a finite number 
of square-integrable CS states. In these integrals, the CS functions, as a function of x a , 
decay exponentially for large x a . Hence the domain of integration is confined to Q a , 
where x a is either finite, or x a « y a as y a — > oo. In this region, as Merkuriev showed 
[2], takes a simple form; it coincides with the channel Coulomb Green's operator 

G a a 3 = G a , (41) 

where G a (z) — (z — Ha)" 1 , and 

H a = H° + v c a . (42) 



Therefore, in calculating the matrix elements in Eq. (39), can be replaced by G a 



Similarly, in calculating (40), V^ s can be replaced by 



U a = vf+v°. (43) 
Consequently, Eq. (J38j) becomes 

{OSY 1 = (GJ- 1 - C/ a , (44) 



where 



and 



G a = a {nvl\; b x by\G a \n'v'l'\'\ b x b y ) a (45) 



U a = a (ni>l\] b x b y \U a \n'i/l'\'; b x b y ) a . (46) 
The \J__a matrix elements can again be evaluated numerically. 

3.4- Matrix elements of G a 

The most crucial point in this procedure is the calculation of the matrix elements G a . 
In our Jacobi coordinates, the three-particle free Hamiltonian can be written as a sum 
of two-particle free Hamiltonians 

= /£, + /£,- (47) 



Thus the Hamiltonian H a of Eq. (42) appears as a sum of two two-body Hamiltonians 
acting on different coordinates 

H a = h Xa + h yaJ (48) 

where h Xa = h Xa + v^(x a ) and h Va = h ya , which, of course, commute. As a result, G a 
is a resolvent of the sum of two commuting Hamiltonians h Xa and h Va . 



Faddeev-Merkuriev integral equations for resonances 



10 



According to the Dunford- Taylor functional calculus, a function of a self-adjoint 
operator h can be defined by 

f{h) = ^-l dz> f{z>){z> -h)~\ (49) 
2m J c 

where C encircles the spectrum of h in counterclockwise direction and / is analytic on 
the area encircled by C. In that way, G a , as a function of the self-adjoint operator h Xa , 
can be written as 

G a (z) = (z-h ya -h Xa Y l 



1 

2tH 



dz (z - h Va - z) 1 (V - h x 

c 



— I dz' g ya (z-z') g Xa {z') , (50) 
J c 



c 

where g Xa {z) = (z — hx^ 1 and g Va (z) = {z — h ya )~ l . The contour C should be taken in 
a counterclockwise direction around the singularities of g Xa such that g Va is analytic on 
the domain encircled by C. Accordingly, to calculate the matrix elements G a , we need 
to calculate a contour integral of the two-body Green's matrices g and g . 

— Vol — £a 

In our case, g Xa is a Coulomb Green's operator with a branch-cut on the [0, oo) 
interval and accumulation of infinitely many bound states at zero energy, while g Va is 
a free Green's operator with branch-cut on the [0, oo) interval. In time-independent 
scattering theory, G a (E) should be understood as G a {E) = lim e ^ G a (E + ie), with 
e > 0. To calculate resonances, we need to continue analytically to e < 0. In this paper, 
we limit our study to energies below the three-body breakup threshold, so dt(E) < 0. 



To examine the analytic structure of the integrand in Eq. (50) let us take e > 0. 
By doing so, the singularities of g Xa and g Va become well separated. Now the spectrum 
of g Xa can easily be encircled so that the singularities of g Va lie outside the encircled 
domain (Fig. [3]). However, this would not be the case for e < 0. Therefore the contour 
C is deformed analytically in such a way that it shrinks to a few lowest bound states and 
the contour opens up and continues along an imaginary line (Fig. [4]). Now, even in the 
e < case (Fig. [5]), the contour avoids the singularities of g Va . Thus, the mathematical 



conditions for the contour integral representation of G a in Eq. (50) are met also for 
resonant-state energies. 

3. 5. The Coulomb-Sturmian matrix elements of the Coulomb Green's operator 

In our system, g Va is a free Green's operator and g Xa is a Coulomb Green's operator. 
Their CS matrix elements can be calculated analytically [TD]. The two-body Coulomb 
Green's operator is the resolvent of the Coulomb Hamiltonian 

gf = {z-hf)-\ (51) 



where 

hP = - . , 

2m V dr 2 r 2 / r 



h? = -¥;($-,- '-^A) + " . (52) 
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8>0 
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Figure 3. The analytic structure of g Va (E + is — z')g Xa (z') as a function of z' , e > 0. 
The operator g Xa (z') has a branch-cut on the [0, oo) interval and accumulation of 
infinitely many bound states at zero energy, while g Va (E + ie — z') has a branch-cut 
on the (— oo, E + ie] interval. The contour C encircles the spectrum of g Xa and avoids 
the singularities of g Va . 



8>0 

* © © © 



t 



Figure 4. The contour of Fig. [3] is deformed analytically such that it shrinks to the 
low-lying bound-state poles of g Xa and the other part is taken along an imaginary 
direction. 
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8<0 
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Figure 5. Even in the e < case, which is needed to calculate resonances, the 
singularities remain separated. Some low-lying poles of g Xa submerge onto the second 
Riemann sheet of g Va , and they are denoted by dotted contour. 



m is the reduced mass, and Z is the strength of the Coulomb potential. In the CS basis 
the operator J = z — hf has an infinite symmetric tridiagonal (Jacobi) matrix structure, 
i.e. all elements are zero, except for the diagonals and off-diagonals, 

jC = 2(z + l + l)(k 2 -b 2 )^- b -Z , (53) 

J^ 1 = -[i(i + 2l + l)]W(k> + #)^ (54) 

and 

J% + i = -[(< + 1)(< + 2/ + 2)] 1 '*{k* + b 2 )^~ b , (55) 

where k = (2mz/h 2 ) 1 ^ 2 . Then, as it has been shown in Ref. [TD], the N x N matrix 
elements of gf are given by 

g C(N) = { jC _ 5m {J C n+i) 2 Cn+i] -1 j (56) 

where J c is the N x N upper left corner of the Jacobi matrix and 



Cn+1 — — — —2 



(b - ik) z (N + / + 2 + i 7 ) 
2 F 1 | -I + i 7 , + 2; N + I + 3 + i 7 ; 



b + ik^ 2 " 



b — ik 



2 F 1 [-l + ij,N+l-N + l + 2 + i 7 ; 



b + ik^ 2 ^ 



b-ik 



(57) 



with 2-^1 being the hypergeometric function and 7 = Z/(mh 2 k). This ratio of two 2-F1 
functions, where the second index in the numerator and in the denominator differ by 
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one, can be represented by a continued fraction [TT], which is easily computable and 
convergent on the whole complex k plane. 

3. 6. Numerical realization of the method 

In this approach for solving the Faddeev-Merkuriev integral equation, the only 
approximation is the replacement of the potentials Va and U a by their respective 
separable forms. We found that good results are achieved when we use N up to 25 — 30 
in the separable expansion for each angular momentum channel. To calculate the matrix 
elements between CS functions, which are, in fact, exponential functions multiplied by 
polynomials, we use Gaussian integration; about 120 — 150 points provide the sufficient 
accuracy. 

The calculation of G a is very accurate. It should be noted first that this 
representation of gj is exact, and its numerical realization, including the evaluation 
of the ratio of two 2 F\ functions by a continued fraction is precise to machine accuracy. 
The contour integral around the poles of g Xa is a projection onto the corresponding 
bound state 

\4>i){4>i\ = ^— < f dzg x A z ) , (58) 
*™ J Ci 

where <pi is the eigenstate belonging to the eigenvalue E iy and Cj is a contour around 
Ei. In fact, the states \<pi) are hydrogenic bound states. We calculated the overlap 



{nl;b\(f)i} using (58) and compared it with the exact result in Maple. We found a 
perfect agreement. We also found that the main contribution to G a is due to the bound 
state poles. The contour integral along the imaginary line behaves asymptotically like 
1/(1 + z') 2 . We adopted the Gauss-rational integration method, and found that about 
50 — 60 integration points provide a sufficient accuracy. This is a significant improvement 
over previous methods which employed as many as 250 integration points in Refs. [U [7] 
to achieve a comparable level of accuracy. 

In order to find those complex zeros of D(E), which are close to the real energy 
line, we have developed the following procedure. We consider an interval along the real 
energy line between two thresholds. Since D(E) is an analytic function of the energy we 
can approximate it with Chebyshev polynomials. We use about n = 12 — 15 Chebyshev 
polynomials. The length of the interval should be small enough that D(E) does not 
change too much and thus the Chebyshev approximation is reliable. The zeros of the 
Chebyshev approximated function is determined by using the eigenvalue method of Ref. 
|12j . Then the rank of the Chebyshev approximation is lowered by one, and the zeros 
are located again. If a zero is a true zero of D(E), the zeros of the rank n and rank 
n — 1 Chebyshev polynomials are close. A similar concept was adopted in Ref. [T5] 
using Pade approximation instead of Chebyshev. We then look for the zeros of D(E) in 
the neighborhood of the zeros of the Chebyshev approximation. We pick three complex 
points, zi, z 2 and z 3 . The location zq of the complex root is estimated by [Ej 
_ zi(z 2 - z^/Djzi) + z 2 (z 3 - z 1 )/D(z 2 ) + z 3 {z x - z 2 )/D(z 3 ) 
Z °~ (z 2 -z 3 )/D(z 1 ) + (z 3 -z 1 )/D(z 2 ) + (z 1 -z 2 )/D(z 3 ) ' 1 ' 
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Table 1. Angular momentum channels I — A used for L — 1 states with natural and 
unnatural parity. The superscript stands for the parity. 



L=l~ 


L = 1 + 


0-1 


1-1 


1-0 


2-2 


1-2 


3-3 


2-1 


4-4 


2-3 


5-5 


3-2 


6-6 


3-4 


7-7 


4-3 


8-8 



Then we make a replacement z — > Zi, Zi — > z 2 and £2 — > £3, and repeat until |z — Zi| < e 
with some small e. If the initial estimation for the zero is good, this procedure converges 
very fast. After some experience, we found this method quite fast and reliable. 

4. Results 

We calculated the resonances of the electron-positronium, e — Ps or e~ — e~ — e + , 
three-body system. Here the two electrons are identical particles, allowing us to use 



the one-component version of the homogeneous Faddeev-Merkuriev equations (21). We 
use atomic units throughout. For the parameters of the cut-off function Q we adopted 
Xo = 30, yo = 35 and v = 2.1. This choice of parameters guarantees that in the energy 
region up to the fifth threshold, there are no spurious solutions. 

The parity of the states is given by P = (—) l+x . If P = (— ) L , the state has natural 
parity, if P = (— ) , the state has unnatural parity. The wave function should be 
antisymmetric with respect to the exchange of the two electrons. If the spin of the two 
electrons couple to S = to form a singlet state, then the wave function is antisymmetric 
with respect to the exchange of electron-spin coordinates, and the spacial part should 
be symmetric. Similarly, if the two electrons couple to S = 1 forming a triplet state, the 
wave function is symmetric with respect to exchange of the spin coordinates, and the 



spatial part of the wave function is antisymmetric. Consequently, in Eq. (21), if S = 
then p = 1 and if S = 1 then p = — 1. 

We present results for total angular momentum L = 1. The angular momentum 
quantum numbers I and A are selected such that I + \ = L . Table [l] shows the angular 
momentum channels used in these calculations. 

In this method, we represent operators on the CS basis, which has one parameter, 
the parameter b. To be economic, we need to find an optimal b, and then we need 
to increase the basis size iV to observe convergence. We found that the results are 
insensitive to varying b over a rather broad interval around b = 0.25. We used b = 0.25 
throughout. Table [2] shows a typical convergence of a resonant-state energy with 
increasing N. From results like this, we can safely infer about three significant digits 
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Table 2. The convergence of a resonant-state energy with increasing N. The results 
are given in atomic units. E r is the real part, Ei is the imaginary part of the energy. 



iV 


E 


Ei 


25 


-0.028958565292 


-0.000000304662 


26 


-0.028959758422 


-0.000000303622 


27 


-0.028960583605 


-0.000000302855 


28 


-0.028961155544 


-0.000000302318 


29 


-0.028961553239 


-0.000000301934 


30 


-0.028961831103 


-0.000000301642 



Table 3. L = 1 natural parity singlet (S=0) and triplet (S=l) resonances in the 
e — Ps system. The thresholds are indicated by empty lines. The resonance energies 
E = E r — ir/2 are given in atomic units. 



5 = 






5=1 




-0.2340 


-0. 


0050i 


-0.248 


-0.003i 


-0.06257 


-0. 


0000000007i 


-0.24 


-0.007i 








-0.23456 


-0.00263i 


-0.0619 


-0. 


0005i 


-0.0619 


-0.0005i 


-0.059948 


-0. 


.0001461i 


-0.0611 


-0.0009i 


-0.0306834 


-0. 


00006229i 


-0.02926 


-0.0000261i 


-0.02896 


-0. 


00000030i 


-0.0281 


-0.000008i 








-0.02794 


-0.00000006i 


-0.018479 


-0. 


00001267i 


-0.0276 


-O.OOOli 


-0.01635 


-0. 


000006i 


-0.01992 


-0.000061i 


-0.0160 


-0. 


0000004i 


-0.0173 


-0.000056i 


-0.01569 


-0. 


OOOOOOli 


-0.0164 


-0.0000007i 


-0.01224 


-0. 


0000063i 


-0.0156 


-O.OOOOOOli 


-0.0107 


-0. 


OOOOli 


-0.0107 


-O.OOOOOli 



for the real part of energy and one or two significant digits for the imaginary part of the 
energy Tables [3] and [4] show the results of our calculations. 

5. Summary 

In this paper, we outlined a solution method for the homogeneous Faddeev-Merkuriev 
integral equations to calculate resonances in atomic three-body systems. We 
approximated the potential terms in the three-body Hilbert space by a separable 
form. This approximation casts the integral equations into a matrix equation and 
the resonances are sought as complex-energy roots of the Fredholm determinant. The 
matrix elements of the three-body channel Coulomb Green's operator were evaluated 
as a complex contour integral of the two-body Coulomb Green's matrices. The use of 
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Table 4. L = 1 unnatural parity singlet and triplet resonances in the e — Ps system. 
The thresholds are indicated by empty lines. The resonance energies E = E r — iL/2 
are given in atomic units. 



s = 


= 






s = 


= 1 






-0. 


.0610 


-0 


.00158i 


-0. 


.061901 


-0. 


0000993i 


-0. 


.061 


-0 


.003i 


-0. 


.060 


-0. 


,002i 


-0. 


.02816 


-0 


,000000002i 


-0. 


.031560 


-0. 


.0000866i 


-0. 


.0277 


-0 


,00003i 


-0. 


.0277 


-0. 


,00006i 


-0. 


.027 


-0 


,0005i 


-0. 


.027 


-0. 


,0005i 


-0. 


.0165 


-0 


,00000005i 


-0. 


.018857 


-0. 


00000690i 


-0. 


.0157 


-0 


,00000003i 


-0. 


.016 


-0. 


000002i 


-0. 


.0108 


-0 


.0000002i 


-0. 


.0125 


-0. 


OOOOli 


-0. 


.0100 


-0 


.0000003i 


-0. 


.0108 


-0. 


00007i 



the Coulomb-Sturmian basis allows analytic evaluation of these matrix elements. We 
found that the contour introduced here is more advantageous than those used in our 
previous publications |H [7j. The method is quite efficient. To achieve good accuracy we 
do not need too many terms in the expansion, only N = 30 in each angular momentum 
channel, and consequently the size of the matrix is relatively small. We performed all of 
our calculations with Mac PC's. We calculated resonances of the e — Ps atomic three- 
body system for total angular momentum L = 1 with natural and unnatural parity. 
We do not believe that there is an ultimate method to calculate resonances. However, 
our results allow us to believe that this solution of the homogeneous Faddeev-Merkuriev 
equations is an accurate and reliable method for calculating resonances in atomic three- 
body systems. 
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